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ABSTRACT 

We propose a solution to the problem of astrometric and photometric cali- 
bration of coronagraphic images with a simple optical device which, in theory, 
is easy to use. Our design uses the Fraunhofer approximation of Fourier optics. 
Placing a periodic grid of wires (we use a square grid) with known width and 
spacing in a pupil plane in front of the occulting coronagraphic focal plane mask 
produces fiducial images of the obscured star at known locations relative to the 
star. We also derive the intensity of these fiducial images in the coronagraphic 
image. These calibrator images can be used for precise relative astrometry, to 
establish companionship of other objects in the field of view through measure- 
ment of common proper motion or common parallax, to determine orbits, and 
to observe disk structure around the star quantitatively. The calibrator spots 
also have known brightness, selectable by the coronagraph designer, permitting 
accurate relative photometry in the coronagraphic image. This technique, which 
enables precision exoplanetary science, is relevant to future coronagraphic instru- 
ments, and is particularly useful for "extreme" adaptive optics and space-based 
coronagraphy. 
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1. Introduction 

Stellar coronagraphs suppress a star's point spread function (PSF) in order to enable 
the imaging of faint objects or structure, such as planets or proto-planetary disks, in orbit 
about the star (Lyot 1939; Sivaramakrishnan et al. 2001). Accurate relative astrometry 
between the star and its apparent companions or extended disk structure is essential to 
establish physical association through confirmation of common proper motion or common 
parallax (Oppenheimer et al. 2001), to determine orbital parameters, and to observe disk 
structure quantitatively. Relative photometry between the star and these other objects 
in the field of view is also necessary to conduct studies of disk or companion physics and 
chemistry. However, a perfect coronagraph will completely remove the core of the star's PSF, 
meaning that measuring its intensity or location in the scientific data is impossible. Lloyd 
& Sivaramakrishnan (2005) show that decentering the on-axis star results in complex PSF 
morphology — tilt errors can even introduce artifacts that look like point sources within a 
few resolution elements of the center of the occulting focal plane mask of the coronagraph, 
even in simple theoretical simulations. We have also observed these effects in real data from 
the optimized, diffraction-limited Lyot Project coronagraph (Sivaramakrishnan et al. 2001; 
Oppenheimer et al. 2003a,b, 2004; Makidon et al. 2005), which is deployed at the Air Force 
AEOS 3.6 m telescope (Roberts & Neyman 2002). This coronagraph is the first to operate 
in the regime of "extreme adaptive optics" (ExAO), with Strehl ratios in the 70-85% range, 
and it has played a role in developing advanced techniques of exoplanetary science. (The 
NICMOS coronagraph on the Hubble Space Telescope is not optimized for coronagraphy 
- the function of its focal plane occulter is to reduce the total flux on the detector, since 
there is no coronagraphic suppression of the stellar background away from the PSF core.) 
It brings early results to a burgeoning field of coronagraphic imaging of nearby stars. In 
Digby et al. (2006) we discuss some of the problems of coronagraphic astrometry uncovered 
when analyzing real data from the Lyot Project coronagraph, including the fact that such 
data is often far more complex than simple theory might predict (Lloyd & Sivaramakrishnan 
2005; Sivaramakrishnan et al. 2005). Over exposures a few seconds to several minutes long 
objects are smeared out on the detector because of slow drifts in the AO system, differential 
refraction effects and field rotation. In addition, the dynamic range of the detector prevents 
simultaneous non-coronagraphic imaging of faint companions with the AO target star. The 
image is further complicated by the motion of speckles, which occurs on all timescales present 
in the obervation. This paper presents a solution to some of the problems discussed in Digby 
et al. (2006). 

We explain how to create controlled fiducial spots in the coronagraphic image which 
pinpoint the location of the star, in addition to providing calibration for accurate relative 
photometry, thus solving the fundamental problem of coronagraphic image calibration. By 
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inserting a reticulate grid of wires into a pupil ahead of the focal plane occulting mask in the 
coronagraph we create predictable images of the central star at desired locations and chosen 
brightness in the coronagraphic image plane. We extend the analytical theory of segmented 
aperture telescope point-spread functions (PSFs) for direct imaging (Chanan & Troy 1999; 
Sivaramakrishnan & Hodge 2001; Yaitskova & Dohlen 2002) as well as that of coronagraphic 
imaging (Sivaramakrishnan & Hodge 2001; Sivaramakrishnan & Yaitskova 2005) to describe 
the effects of our proposed reticulate grid of wires in the coronagraphic case. 

Our methods, though relevant to any coronagraphic observation, are particularly ap- 
plicable to space-based high dynamic range coronagraphy dedicated to the extremely chal- 
lenging task of finding and characterizing Earth-like exoplanets, as our approach does not 
disturb the extremely smooth, well-corrected wavefront that such searches demand. Our 
technique does not require exquisite calibration of deformable mirrors, for instance, nor does 
it require any change in a deformable mirror shape before, during, or after an observation. 

We note a simultaneous but independent study by Marois et al. (2006), where a pupil 
plane before the focal plane mask is also modified in order to create images suitable for 
photometric and astrometric purposes. 

§2 elucidates the theoretical basis for our calculations of telescope PSFs. §3 describes 
the perfect coronagraph used to illustrate the utility of this new theoretical construct, as 
well as the difficulties in proper numerical simulations of the reticulate pupil mask. In §4, 
we extend the theory to a more practical level, addressing the effects of finite bandwidth in 
astronomical imaging, imprecise positioning of optics, residual uncorrected wavefront errors, 
and atmospheric differential refraction. §5 describes the practical application of the theory, 
including some details on how to design reticulate pupil masks, and how to calibrate and 
reduce coronagraphic data to derive precision relative astrometry and photometry. 

2. Point-spread function theory 

In this section we describe our monochromatic Fourier optics formalism. We recollect 
that a plane monochromatic wave traveling in the z direction in a homogenous medium 
without loss of energy can be characterized by a complex amplitude E representing the 
transverse (e.g., electric) field strength of the wave. The full spatio-temporal expression for 
the field strength is Ee 1 ^^^, where lu/k = c, the speed of the wave. We do not use the 
term field to denote image planes — the traditional optics usage — we always use the term 
to denote electromagnetic fields or scalar simplifications of them. The wavelength of the 
wave is A = 2ir/n. The time-averaged intensity of a wave at a point is proportional to EE*, 
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where the average is taken over one period, T = 2n/u, of the harmonic wave. The phase 
of the complex number E represents a phase difference from the reference phase associated 
with the wave. A real, positive E corresponds to an electric field oscillating in phase with our 
reference wave. A purely imaginary positive value of E indicates that the electric field lags 
by a quarter cycle from the reference traveling wave. Transmission through passive, linear 
filters such as apertures, stops and apodizers is represented by multiplying the field strength 
by the transmission of these objects which modify the wave. Again, such multiplicative 
modification is accomplished using complex numbers to represent phase changes forced on 
the wave incident on such objects. 

We assume that the Fraunhofer approximation describes our imaging system. Thus 
image field strengths are given by the Fourier transform of aperture (or pupil — we use the 
two terms interchangeably) illumination functions, and vice versa. 



Gaps between segments are analogous to the wires in our grid, and our wire grid spacing 
replaces the segmentation periodicity. We use the notation of segmented apertures in order 
to facilitate comparison of this work with existing work on segmented aperture PSFs, but 
refer to wire thickness and spacing rather than segment gap and width. We re-derive the form 
of square-segmented aperture PSFs in both direct and coronagraphic imaging analytically, 
extending the approach of Sivaramakrishnan & Hodge (2001); Sivaramakrishnan & Yaitskova 



The full aperture of the telescope is described by a function A(x), which, for unapodized 
apertures, is unity over the aperture and zero elsewhere. The position vector x = (x, y) 
denotes location in the aperture plane, in units of the wavelength of the light. We develop 
the expression for the monochromatic PSF with a perfect wavefront first, and then calculate 
the coronagraphic PSF, assuming that the scalar wave approximation and Fourier optics 
provide an adequate description of PSF formation. We treat aberrated wavefronts in Section 



The aperture illumination function — the complex amplitude in the aperture in response 
to a unit (field rather than power) strength incoming wave — for the reticulated aperture, 
can be written as 



2.1. The reticulated aperture point-spread function 



(2005). 



4.3. 



4,(x) 



= A(x) x 

(l-lll((x-x )/d) / d*U(x/g)) x 
(l-lll((y-y )/d) / d*U(y/g)). 



(1) 
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Here d is the wire grid's spatial periodicity, g the wire width; * denotes the convolution 
operation. We make use of Fourier analytic notation, techniques and results which can be 
found in e.g., Bracewell (1986). The one- dimensional shah or comb function III is defined 
by 

oo 

11%) = S(u-n), (2) 

n=— oo 

where 5(u) is the Dirac 5-function, and 

U(x) = 1 for \x\ < 1/2, 
Il(x) = elsewhere. 

Normalization by 1/d is required to ensure the correct impulse strength in each 5-function in 
equation (1). We offset the origin of the grid from that of the full aperture by an arbitrary 
vector shift of (x ,y ), restricting x Q and y Q to be less than the grid size d (without loss of 
generality). This lateral offset of the reticulate pattern is a detail that segmented aperture 
telescope PSF calculations (e.g., Chanan & Troy (1999); Sivaramakrishnan & Hodge (2001); 
Yaitskova et al. (2003); Sivaramakrishnan & Yaitskova (2005)) do not need to consider, but 
is relevant to our case, as it enables analytical tolerancing of the in-plane mechanical stability 
and placement of the wire grid. 

The complex amplitude of the image resulting from the aperture illumination function 
A s is its Fourier transform, 

a s (k) = a(k) * 

(S(k x ) - g lll(dk x ) e tk * x ° smc(gk x )/d) * 

(5(k y ) - g lll(dk y ) e lkyy ° smc(gk y )/d), (4) 

where k = (/ci,/^) is the image plane coordinate in radians. We adopt the convention of 
changing the case of a function to indicate its Fourier transform. We refer to this complex- 
valued field a as the "amplitude-spread function" (ASF), by analogy with the PSF of an 
optical system. The ASF is the electric field in the image plane. The PSF is the square of 
the absolute value of the ASF (equation (7), below). Using the definitions 

oo 

G x (k) = - e ik * x ° a{k x -n/d,k y ) smc(gk x ) 

n=—oo 
oo 

G y (k) = -e ih « y ° ^ a(k x ,k y - m/d) smc(gk y ) 

m=— oo 

oo 

G xy (k) = e i{k * x ° +k y yo) a(k x -n/d,k y -m/d) 

n,m=~oo 

x sinc(gk x ) smc(gk y ), (5) 
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the expression for the ASF of the reticulated pupil is 

a s (k) = a(k)+e(G x (k) + G y (k)) + e 2 G xy (k) 

(6) 

We use the dimensionless number e = g/d to characterize the wire width relative to the wire 
spacing. The corresponding PSF is 

Ps (k) = | a |2 + e 2(| Gx |2 + | Gi/ |2 )+e 4| Ga:w |2 

+ e (a(G x + G* y )+a*(G x + G y ))+0(e 3 ), (7) 

We explicitly write the e 4 \G xy \ 2 term in spite of the presence of a 'catch-all' 3 rd order re- 
mainder term 0(e 3 ) because the former has peaks that could on occasion be significant 
in the coronagraphic PSF. At the origin the expression for a s in equation (6) reduces to 
(1 — e) 2 a(0, 0), so the peak value of the reticulated aperture PSF is (1 — e) 4 p(k), or 
~ (1 — 4e) p(k) if e 1 (by the binomial approximation) where p = aa* is the full aper- 
ture PSF. The accuracy of this approximation affects our ability to estimate photometric 
calibration analytically, but in any case the satellite PSFs will require their own calibration 
methods. In practical applications e will likely be far smaller than the value we use in our 
numerical example below (see §5.1), so use of the binomial approximation produces useful 
analytical insight by enabling a simple theoretically predicted expression for the value of the 
PSF maximum. 



2.2. Components of the point-spread function 

The PSF in equation (7) can be decomposed into a few physically significant compo- 
nents. The PSF has a bright central core of (1 — 4e) p(k). By inspection of the e 2 \G x \ 2 and 
e 2 \G y \ 2 terms we see that the PSF possesses 'satellite' copies of p arranged along the arms 
of a cross, on two linear arrays at a spacing X/d on the k x and k y axes, with an overall taper 
by sinc(gk x ) and sinc(gk y ) respectively. (Fig. 1, panels 1 and 3 are monochromatic, panel 
5 shows a 20% bandpass polychromatic PSF). The e 4 (G^l 2 term in equation (7) yields a 
sea of copies of e 4 p in a two-dimensional grid off the axes in k-space, also the same X/d 
spacing. This 0(e A ) component is not visible in the direct images in Fig. 1 (panels 1, 3, and 
5) because they are too faint to be seen with the grey scale used in these panels. However, 
they are apparent in the coronagraphic PSFs using the same reticulated aperture (Fig. 1, 
panels 2, 4, and 6), where the central core has been suppressed. There are first and higher 
order "fringing effects" in the PSF (Fig. 3). These are described by cross terms in equation 
(7). The first order fringing term is e(a(G* x + G*) + a*(G x + G y )). When x Q or y Q are non-zero 
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these fringing effects can introduce a complicated (but still centro-symmetric, since A s (x) is 
a real function) morphology in the image. The phenomenon of speckle pinning (Bloemhof 
et al. 2001; Sivaramakrishnan et al. 2002; Perrin et al. 2003; Aime & Soummer 2004) will 
cause complex interactions with residual phase aberrations in the brighter parts of the im- 
age, so where these fringes overlap with the bright satellite peaks their effects on astrometry 
and photometry will be more severe than in the darker areas inbetween satellite peaks. We 
discuss this further in Section §4.3. 

If there are many grid spacings across the full aperture, X/D <C X/d, and the fringing 
between various components of the ASF is greatly reduced in the PSF. In this case the 
PSF in equation (4) separates out into obvious physical components deriving from the full 
aperture ASF and the squares of the absolute values of G X) G y and G xy . This entire pattern 
of full aperture PSFs is weighted by a broad 2-dimensional sine function envelope with a 
characteristic scale X/g, which is the widest angular scale in the ASF. This scale derives 
from the inverse of the smallest physical distance in the reticulate aperture, viz., the wire 
thickness. This overall weighting is almost constant at angular distances of ~ 0.1 X/g from 
the center of the PSF, simply because of the properties of the sine function. Thus both 
sets of satellite PSFs (on the k x , k y well as off them) have more or less uniform 

intensities near the core. This morphology is visible in Fig. 1. As e — > 0, p s (k) reduces to a 
single central PSF p(k) at the origin — satellite spots and fringing disappear completely, as 
expected, when the wires vanish. 

Alternate geometries can be created with different spacings in the x and y axes when 
using specialized entrance apertures, such as that currently envisaged for the Terrestrial 
Planet Finder. In point of fact any periodic two-dimensional wire grid can be used. 

Before we discuss the utility of forming PSFs of this nature, we must derive the corona- 
graphic PSF. This is done in the next section. §4, then, deals with the effects of finite spectral 
bandwidth, an imprecisely placed reticulate pupil mask and residual phase aberrations in 
the entrance pupil. 

3. Coronagraphic PSF on a reticulated aperture 

A telescope aperture is described by a transmission function pattern -A(x), where x = 
(xi,x 2 ) is the location in the aperture, in units of the wavelength of the light (see Fig. 2). 
The corresponding aperture illumination describing the electric field strength in the pupil (in 
response to an unaberrated, unit field strength, monochromatic incident wave) is E\ = A(x). 
The field strength in the image plane, E B = a(k), is the Fourier transform of E A where 
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once more, k = (k±, k 2 ) is the image plane coordinate in radians. Because of the Fourier 
relationship between pupil and image fields, k is also a spatial frequency vector for a given 
wavelength of light. The PSF is p = aa*, recalling our convention of changing the case of a 
function to indicate its Fourier transform. We, again, refer to a as the ASF. We multiply the 
image field Eb by a mask function m(k) to model the focal plane mask of a coronagraph. 
The image field immediately after this mask is Eq = m(k)E B . The electric field in the 
re-imaged pupil after the focal plane mask (the Lyot pupil) is E D , which is the Fourier 
transform of E c . We use the fact that the transform of the image plane field E B is just the 
aperture illumination function Ea itself, so the Lyot pupil field is E D = M(x)*£U, where * 
is the convolution operator. 

If the Lyot pupil stop transmission is iV(x), the electric field after the Lyot stop is 
E E = N(x)E D . The transform of this last expression is the final coronagraphic image 
field strength: E F = n(k)*[ra(k) E B \- Sivaramakrishnan & Yaitskova (2005) describe the 
structure of the field strength Ed in the Lyot plane located at D. We use their results 
to analyze the final image plane Ep for an ideal ("perfect") coronagraph on a reticulated 
circular aperture, in order to use its features as fiducial locations in the final image for the 
purpose of high precision astrometry and photometry. 

3.1. The ideal coronagraph 

The band-limited coronagraph (BLC) design is mathematically "perfect" in that simple 
Fourier optics modelling predicts that it will prevent all incoming, on-axis light from reaching 
the final coronagraphic focal plane (Kuchner & Traub 2002). Its mathematical simplicity 
makes it useful for the purpose of elucidating the way Lyot-style coronagraphic PSFs are 
affected by various factors, such as tip-tilt (Lloyd & Sivaramakrishnan 2005), higher order 
phase aberrations (Sivaramakrishnan et al. 2005), the presence of secondary mirror support 
vanes or "spiders" in the pupil (Sivaramakrishnan & Lloyd 2005), the secondary obstruction 
itself (Sivaramakrishnan et al. 2004), or inter-segment gaps (Sivaramakrishnan & Yaitskova 
2005). The BLC design also enables an analytical treatment of various effects, which can 
then be studied numerically for particular instruments, or on related designs such as the 
Apodized Pupil Lyot Coronagraph (Soummer et al. 2003; Soummer 2005; Sivaramakrishnan 
& Lloyd 2005). 

The purpose of this exercise is that although the BLC may not be used in any particular 
coronagraphic instrument, the theory allows us to examine the behavior of the reticulated 
mask independent of the particular choice of coronagraphic starlight suppression. In this 
way, our analysis is not complicated by the difficulties of starlight suppression, by the simple 
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fact that we assume it is perfect in this analysis. 

The BLC design is easily understood if we introduce a "mask shape function" w with 
the definition 

m(k) = 1 -io(k), (8) 

which means that M(x) = 5(x) — W(x). We use this in the expression for the Lyot plane 
field strength Ed, noting that since the focal plane mask is opaque at its center, w is unity 
at the origin, so W is a function whose two-dimensional integral over the whole of the pupil 
plane is unity Thus we can write the electric field at the Lyot plane as 

E D (x) = A*(8(x) - W(x)) 
= A - A*W(x). 

If w is a band-limited function with bandpass b, then there is a minimum positive value of 
b such that the mask function's Fourier transform, W, satisfies the property 

W(x) = if |x| > b. (10) 

Thus in the BLC design W (the Fourier transform of the mask shape function) has compact 
support in its domain (the Lyot pupil plane). 

In the case of a completely unaberrated coronagraphic optical train, the above-mentioned 
properties of w (or W) result in the Lyot pupil field at any point further than 2b from a pupil 
edge being identically zero. We must also invoke the fact that the area under the convolution 
of two functions is the product of the areas under the two constitutent functions to see that 
this is indeed so. 

In our numerical examples of coronagraphic images, which are all unobscured circular 
aperture systems, we use a BLC with a mask profile 

m(k) = l-jinc 2 (|k|/A; ), (11) 

where jinc(a;) = 2J 1 (x)/x. Here J\{x) is the Bessel function of the first kind, with index 
1. The constant k Q is chosen so that the first zero of the jinc function lies 12 resolution 
elements (A/-D) from the center of the opaque focal plane mask. Since this is about 10 times 
as wide as the direct image's PSF, the bandpass of this jinc 2 mask is approximately D/5 
(the bandpass of a product of two functions being the sum of the bandpasses of the two 
functions, the bandpass of our jinc function is D/10). Therefore the optimally-sized Lyot 
stop for this "fourth-order" focal plane mask is 3D/5 in diameter (the order here refers to 
the local angular intensity profile of the suppressed or nulled area of the coronagraph) . 



(9) 
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3.2. Numerical simulations 

Our numerical example in Fig. 1 uses a 128-sample wide binary pupil in an array of 1024 
elements on a side, wires 2 samples wide, with 8 wire gaps across the pupil in the x\ and X2 
directions. Those in Figs. 3, and 4 use a 512-sample pupil embedded in a 2048 x 2048 array, 
wires also 2 samples wide, and 32 samples between wires, resulting in 16 wire spacings across 
the pupil. Our fast Fourier transforms in double-precision (64-bit) numbers are implemented 
in the FFTW package (Frigo & Johnson 1997) wrapped in Python's numarray package 
(Greenfield et al. 2002; Greenfield et al. 2003). 

The lattice of satellite spot locations in our image plane are at (lQn\/D, lQm\/ D), 
where n and m are integers (at least one of them being non-zero). Two pupil samples across 
each wire are barely adequate to sample the wires, although this suffices for illustrative 
purposes. However, there is a clear separation of scales between our wire thickness and wire 
spacing with these parameters. The wire thickness to wire spacing ratio, e, is 1/16, so our 
first few satellite spots along the k x and k y axes are e 2 = 3.91 x 10~ 3 (6 magnitudes) fainter 
than the central PSF core in the direct image. Our coarsely-sampled simulations confirm 
the theoretically calculated intensity ratios to within 4%. 

However, the spots in the coronagraphic image are a factor of 5/3 wider (in linear extent) 
than the direct image's spots, because the optimized Lyot stop is 3/5 as large as the entrance 
pupil. This also means that the spots, while the same angular distance from the core, are 3/5 
closer in terms of the effective resolution of the system because of the Lyot stop undersizing 
in this Lyot coronagraph design. This change in size is dependent on the choice of starlight 
suppression, and is not a property of the reticulate pupil mask. Computers available now 
are prohibitively slow at simulating a two-dimensional grid of fine wires across a pupil at the 
scales of interest to the exoplanet hunter, so analytical estimates are useful for understanding 
the behavior of this configuration. For example, to create fiducial calibration spots on the 
order of the brightness of a planet, so that the calibration does not overwhelm the planet's 
brightness, one would want fiducial spots of the order of 15 magnitudes fainter than the star, 
which is approximately the if -band flux ratio between a 40 Myr old Jupiter-mass companion 
and its host star (Burrows et al. 2004). Thus, the wires would need to be one thousandth 
as thick as their spacing. Placing two numerical samples across the wires in this case would 
entail 32768 x 32768 pupil arrays. Typically coronagraphic computations require oversizing 
these arrays by another linear factor of 4 to 8 in order to provide sufficiently fine sampling in 
the pupil plane. Thus numerical FFT methods will need to handle arrays 196,000 samples on 
a side. While direct PSF calculations involving the semi-analytical 'grey pixel approximation' 
do not require such array oversizing (Troy & Chanan 2003), this approximation may not be 
accurate for coronagraphic modelling. 
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The results of our crude calculations of reticulate mask behavior are shown in Fig 1. 
The left hand column (panels 1 and 3) are direct images of a star with the reticulate mask 
in place. The right column (panels 2 and 4) shows coronagraphic images. Panels 1 and 2 
are at the shortest wavelength of a typical astronomical filter's bandwidth (~ 20%), while 
panels 3 and 4 show images made at the longest wavelength. Panels 5 and 6 are broad- 
band images. We now describe the derivation of these broad band images, as well as the 
technique's sensitivity to optic placement and residual wavefront error. 

4. Broad-band theory and tolerancing 

The theory and numerical models we present in the preceding section apply to the case 
of monochromatic light. A comparison of the two upper sets of panels in Fig. 1 shows that 
the fiducial spots to be used for astrometric and photometric calibration are at (angular) 
positions in the focal plane that are a function of wavelength. Indeed, the spectral width 
of an observing bandpass induces radial smearing of the fiducial spots. These spots will be 
miniature, low-resolution spectra of the central, obscured star. Intuitively, this is due to the 
fact that the reticulate mask acts as an extremely coarse grating inserted into the beam. 

For the purposes of astrometry and photometry, we are interested in the four closest 
satellite PSFs located on the k x and k y axes. These four spots lie at (±\/d, 0) and (0, ±\/d) 
relative to the central star. At this point we ignore any higher-order spots, assuming that 
the coronagraph that uses this technique will be designed so that this first set of spots are 
placed as far away from the central star as possible, while still remaining in the field of view, 
so as to obscure as little of the potential search area around the star as possible. In §4.2 we 
argue that there are benefits that accrue if the satellite spots lie within the AO control area 
of the wavefront correction system used in tandem with the coronagraph. 

4.1. Satellite spot intensity 

The four spots closest to the star from a reticulate pupil with N = D/d wire spacings 
across the full pupil (of diameter D) will lie at (±N\/ D,0) and (0,±NX/D) in the image 
plane at any particular wavelength A. If the observing bandpass, with central wavelength 
A c , has a fractional bandwidth f3 = S\/X c , the radial smearing of one of these spots is 
N5X/D = N(5\/\ c )\ c /D, or about Nf3\ c /D. The spots themselves extend over the coor- 
dinates (±N(3\ C /2D, 0) and (0, ±N(3\ C /2D). Approximating the equivalent width of PSFs 
in the entire spectral band by that of the PSF at the central wavelength of the bandpass, 
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A c , we see that each radial smear covers approximately 1 + N(3 square resolution elements. 
Fig. 1 (panels 5 and 6) shows a numerical simulation of a broad-band PSF with a 20% 
bandwidth. We assume a perfect detector and a perfectly flat stellar spectrum, with sharp 
cutoffs at both the long and short wavelength ends of the filter pass band in these figures. 

These spot widths will depend on the type of coronagraph. In classical unapodized 
coronagraphs {e.g., Sivaramakrishnan et al. (2001); Sivaramakrishnan & Yaitskova (2005)) 
the satellite image width is determined by the effective pupil diameter of the system after 
the undersized Lyot stop. This affects the final image size by effectively reducing D, though 
the reticulate mask will act on the entrance pupil of size D. Lyot pupil undersizing affects 
the width of the satellite spots rather than their angular location in the focal plane. This 
detail must be taken into account when estimating the satellite PSF brightness, for example, 
in order to design a reticulate mask for a given application. 

The total flux in the four closest satellite images is 4e 2 of the stellar on-axis flux in the 
absence of a focal plane mask, but through the undersized Lyot stop. Using this for accurate 
photometry is discussed in 5.4. We show in §5.3 that the fact that each smeared spot is a 
stellar spectrum does not affect the precision of the astrometry. 

4.2. Sensitivity to in-plane positioning 

In order for this reticulate mask to be useful for astrometric and photometric measure- 
ments, we must show that slight misalignments of the reticulate mask will not affect whether 
the star can be accurately located, nor will they grossly affect the fiducial spot intensities. 
Here we show that the characteristics of the spots relative to the star are in fact determined 
by the reticulate mask characteristics, not its position in the pupil plane. We quantify the 
limits of this statement below. 



4-2.1. Photometric effects 

The number of wire grid spacings across the aperture, N, determines the relative magni- 
tude of the cross terms caused by interference between satellite peaks. This parameter feeds 
into the PSF structure through d = D/N, thus affecting the separation between satellite 
spots. We have already referred to this interference between satellite peaks as fringing. 

A coronagraph suppresses the central peak of the PSF, but lets the satellite peaks 
through to the coronagraphic focal plane (see (Sivaramakrishnan & Yaitskova 2005) for 
details). Therefore in the coronagraphic case fringing only occurs between satellite peaks. 
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Fringing between nearest neighbor satellite peaks on the k x or k y axes yield field strength 

3 

amplitudes within the range ea(0, 0) ±ea(0, 0)N~2. We have used the fact that at a distance 
of several resolution elements from the origin, a(k) ~ for a circular entrance aperture 
telescope by virtue of the asymptotic behavior of Bessel functions {e.g., Abramowitz & 
Stegun (1970)). 

Thus the absolute PSF intensity variation of our satellite peak due to fringing with its 
nearest neighbor is of the order of e 2 a(0, 0)(l±2iV _ 2). Remembering that the coronagraphic 
satellite peaks on the axes have an intensity e 2 p(k), the fractional change in satellite PSF 

3 

brightness is therefore of the order of 2N~z. 

If iV ~ 25, the resultant fringing is of the order of 1/64, or 1.6%, of the satellite image's 
peak intensity. For N = 50 this reduces to fluctuations of at most 0.6%. Examples of the 
fringe morphology are shown in Fig 3. 

As a result, arbitrary decentering of the reticulate mask affects photometric calibration 
by no more than 1.6% in our example. Although this places a limit on the photometric 
precision attainable with this method, far superior performance can be achieved by requiring 
a modest level of alignment between the reticulate mask and the pupil. In any case, 1.6% 
photometry of any exoplanet is certainly considerably more accurate than any published 
exoplanet photometry to date. 

When phase aberrations are present the ASF is given by a s (k) * FT(e 1 ^). If satellite 
images are separated by more than the size of a seeing disk, their peaks will track the 
Strehl ratio of the core closely. This can be seen by the fact that a s (k), the perfect ASF, 
is basically a set of (^-functions on a grid with spacing \/d. This ASF is convolved by the 
Fourier transform of e l<t> , the latter being the seeing disk. Given that the value of FT(e^) 
at the origin tracks the Strehl ratio exactly, each satellite peak will do so also, as long as 
FT(e t< t > ) drops to a negligible value at a distance \/d from its origin, viz., the seeing disk 
is significantly smaller than the satellite spacing. When satellite images lie in the wings of 
the PSF the situation is more complicated because the closest satellite peaks will lie on a 
pedestal of a seeing halo around the core. Of course the addition of the halo and the wings 
is coherent rather than incoherent. This situation needs further analysis to understand the 
photometric correction of satellite intensities due to the core PSF's halo. Given that an 
ExAO system reduces the parent full aperture PSF's intensity to some noise floor within 
the AO control area, and the satellite peaks must rise above this noise floor in order to be 
identified and measured, the AO-corrected core will not affect the satellite PSF photometry 
by more than the ratio of the satellite intensity to the noise floor induced by residual phase 
errors. The effects of the nearest neighbor satellite peak will be down a factor of e 2 from the 
effect of the core on the satellite. This will have a negligible effect on the photometry of the 
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satellite images. 

4-2.2. Astro-metric effects 

From a simple theoretical standpoint, the locations of the spots will always remain 
symmetric about the occulted star's position on the detector in the absence of phase errors 
and differential refraction through the atmosphere or transmissive optical elements. 

A coronagraphic imaging survey designed to detect exoj up iters over a one to four year 
period requires astrometric precision on the order of a milliarcsecond, or X/50D on an 8 m 
telescope in the if-band. Terrestrial Planet Finder science goals and unpublished studies 
conducted by both the Gemini and VLT planet imaging teams require 1 to 10 mas precision 
in relative astrometry for discrimination between background objects and real companions 
of nearby stars when orbital motion is included. There is no reason why ~ l/100th of a 
PSF width would not be achievable, particularly with four fiducial spots. The astrometry 
of the putative planet relative to the central star will be predominantly determined by the 
statistical significance of the planet detection. 

A tilt of the mask to the beam results in a diminution of the effective value of d, since the 
projected separation in the pupil plane will contain a factor of the cosine of the tilt angle. A 
tilt of 12 degrees is required to produce a projected separation change of a factor of 1 — 1/50, 
so mask tilt in a collimated beam does not need to be tightly controlled. If the reticulate 
mask is allowed to rotate as well as translate, our method is not in principle sensitive to the 
azimuthal orientation of the grid either, as long as our data reduction algorithms allow for 
such rotation of the pattern of fiducial spots. 

Ground-based systems using the reticulated pupil mask for astrometry will also have 
to contend with differential refraction from the atmosphere. Even with an atmospheric 
differential refraction corrector (ADC) there will be residual differential refraction. This will 
distort the small spectral smears of the satellite spots. We have not estimated the magnitude 
of this effect, as it depends on particular details of ADC design, as well as atmospheric 
properties (Roe 2002). By observing a star at different zenith distances the residual errors due 
to atmospheric differential refraction can be calibrated out to some extent. This calibration 
may also be dependent on the spectrum of the star. Furthermore, too much reticulate mask 
rotation is likely to complicate the calibration of atmospheric differential refraction. 

We conducted numerical simulations of in-plane placement errors. The nine corona- 
graphic images in Fig. 3 show the effect of random decenterings of the reticulate mask 
within the unit cell of the grid pattern (a square of side d). We use random placements to 
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avoid potential 'resonances' between predetermined offsets that are themselves in a periodic 
grid within the unit cell of the wire grid. Simulation of many more random placements, 
in addition to systematic incremental grid decentering, indicated that the first nine cases, 
shown here, produce generic behavior. Given the centro-symmetry of the coronagraphic im- 
ages (which ensues since the pupil illumination function is a real function), no change in the 
estimation of the center of the pattern will result if the astrometric data reduction method 
utilizes this fact. In the absence of aberrations and photon noise, in-plane grid placement 
errors make no intrinsic contribution to astrometric errors, unless image distortion calibra- 
tion errors are significant enough to destroy the symmetry inherent in using a reticulate wire 
grid in the pupil. Thus image distortion needs to be characterized to at least the same level 
as the required astrometric accuracy demanded by the science. 

More sophisticated wave propagation and geometric distortion models are needed to 
quantify astrometric sensitivity to lateral placement errors of the grid in any particular 
instrument. In the distortion- free Fraunhofer regime lateral placement of the reticulate grid 
has no effect on the estimate of the occulted star's position in the coronagraphic image. 

4.3. Sensitivity to residual phase aberrations 

As Perrin et al. (2003) point out, when pupil-plane phase aberrations are the dominant 
aberrations, the PSF of an imaging system can be decomposed into an infinite but absolutely 
convergent series involving a combination of the perfect system's response and the Fourier 
transform of residual phase aberrations over the pupil. This work suggests that the fiducial 
spots' morphology can be changed in complicated ways by these interactions, even in the 
case when the pupil-to-image relation is described by simple Fourier optics. In light of these 
facts, we examine how residual phase aberrations might affect the astrometric performance 
of the fiducial spots created by a reticulate grid in the pupil. 

The left column of Fig. 4 (labelled "No WFE" ) shows three realizations of the unaber- 
rated coronagraphic response of a perfect coronagraph, with the grid origin placed at three 
random locations in the entrance pupil of the telescope. This simulates in-plane placement 
errors of the mechanism holding the grid. 

This coronagraphic configuration is illuminated with aberrated but ExAO-corrected 
on-axis starlight. We simulate ExAO correction by high-pass filtering a realization of a Kol- 
mogorov spectrum phase screen, with D/r Q = 15. The high-pass spatial frequency filter 
mimics the action of a perfect ExAO system equipped with spatially-filtered wavefront sens- 
ing (SFWFS) (Poyneer & Macintosh 2004). We generate a phase screen using the method 
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described by Lane et al. (1992), and our numerical filtering procedure follows that of Sivara- 
makrishnan et al. (2001); Poyneer & Macintosh (2004); Makidon et al. (2005). Our spatial 
filter blocks all spatial frequencies below the filter's cutoff, but passes all signal at spatial 
frequencies above the cutoff. The cut-off spatial frequency is determined by the number of 
actuators across the telescope diameter (as long as wavefront sensing is performed on at least 
as fine a spatial sampling as the control). This frequency is referred to as the ExAO control 
frequency (Perrin et al. (2003); Poyneer & Macintosh (2004) discuss this topic further). The 
wires in our grid are spaced by d = D/32. We then simulate two different ExAO systems, 
one with 24 actuators across the entrance pupil (iV ac t = 24), which places the fiducial spots 
well outside the ExAO control area, and another in which N act = 40 to ensure the spots are 
located inside the ExAO control area. This allows us to evaluate how residual speckles in 
the final image affect the utility of the satellite calibrator spots. Coronagraphic PSFs from 
three realizations of grid placement in the pupil of these two AO systems are shown in Fig 
4. 

After the spatial filtering simulating the ExAO, we rescaled the our phase samples to 
ensure that the variance of the filtered phase array is 0.05 radians 2 , to ensure a Strehl 
ratio similar to that expected for real ExAO systems being designed now. The Marechal 
approximation (Marechal 1947) states that this variance is the Strehl ratio reduction, so 
a 95% Strehl ratio image is placed on our coronagraph's focal plane mask. Thus 5% of 
the energy in the AO-corrected wave is distributed in remnant speckles. This energy leaks 
through the coronagraph (Sivaramakrishnan et al. 2001; Lloyd & Sivaramakrishnan 2005; 
Sivaramakrishnan et al. 2005), and affects the shape and symmetry of the fiducial spots in 
the coronagraphic image plane. 

The middle column of Fig. 4 (labelled "iV act 24" ) shows the same three realizations of 
pupil grid geometry as in the "No WFE" column, but with our Kolmogorov-spectrum phase 
screen corrected by an ExAO system with 24 actuators across the entrance pupil. With this 
arrangement the fiducial spots lie in the midst of uncontrolled residual speckle. The peak 
intensity of the four satellite spots in this image is at a level of ~ 3 x 10~ 4 of that of the 
unocculted star. We use a logarithmic grey scale between this peak and 1/251 of this peak 
(6 magnitudes of dynamic range) in the figure. The darkest regions of our images are 14.8 
magnitudes fainter than an unocculted star. At this display stretch we see that the fiducial 
spots lose their simple morphology, as well as their symmetry, when immersed in a sea of 
residual speckles. 

The right column of Fig. 4 (labelled u N act 40") shows the same three realizations of 
pupil grid geometry as in the left column, but with 40 actuators across the entrance pupil. 
Now the fiducial spots lie inside the darker ExAO control area, and are less prone to the 
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deleterious effects of remnant speckles than the A^ act = 24 case shown in the middle column. 
We discuss these practical matters in the next section. 

5. Application of the Theory 

The previous sections describe a new reticulate mask that produces fiducial spots for 
astrometric and photometric calibration of coronagraphic images. Using the Fraunhofer 
approximation of Fourier optics, we have shown analytically, as well as with numerical simu- 
lations, how these spots are formed in both monochromatic and broad-band light. We have 
also demonstrated that their brightness and spacing is easily controlled, and that extremely 
precise positioning of the reticulate mask is not required for the use of the satellite image 
spots for astrometry or photometry, and furthermore, that our technique is usable with real 
ExAO corrected starlight. In this section we explain how to design the mask, how it would be 
used, and how astrometry and photometry can be derived from real data using our reticulate 
pupil mask. 

5.1. Reticulate mask design 

As with many new concepts, there are several instrument-specific details that must 
be addressed by an initial design. Further refinements of any early design will doubtless 
occur during the manufacture and use of a reticulate pupil mask on a telescope. Practical 
problems such as the quality of pupil image at the grid will have to be studied — this 
particular problem appears to be tractable, since coronagraphs typically require good pupil 
images, and beams in coronagraphic instruments have high focal ratios, and, thus, slowly- 
converging beams (Oppenheimer et al. 2003a, 2004; Macintosh et al. 2004; Mouillet et al. 
2004). 

Residual speckle noise in a coronagraphic image sets the contrast ratio floor that the 
instrument can achieve unless post-processing methods (e.g., Marois et al. (2000); Marois 
(2004); Marois et al. (2004)) can suppress this speckle noise further. Here we assume coro- 
nagraphic data is not processed further in such ways. Hinkley et al. (2006) report residual 
speckle noise in the if-band at the edge of the AO control area (17X/D from the PSF core, 
or 1".5 at the if -band central wavelength) in the Lyot Project coronagraph at ~ 11 magni- 
tudes fainter than the central star. We might therefore wish to create fiducial spots that are 
clearly visible, but not overpoweringly bright, near the edge of our usable AO-corrected field 
of view. The location of the fiducial spots is determined by the periodicity of the grid of 
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wires relative to the full pupil diameter, and the spots' brightness is set by the wire thickness 
relative to the wire spacing, e = g/d. With a 100 mm pupil at the beam capture mirror BCM 
of the Lyot Project coronagraph, a wire spacing of 5.88 mm will place the fiducial spots at 
17X/D (l'/5 at if -band center) along the grid's symmetry axes in the image plane. In order 
to make the spots 10 magnitudes fainter than the central star, we require a wire thickness 
of 1/100 of their spacing, viz., 58.8 fim diameter. 

In summary, the relations describing the design are that the four closest and brightest 
satellites spots are placed (D/d)(X/D) = X/d distant from the central star (d being the wire 
grid spacing as projected on to the entrance aperture), and their intensity is a fraction (g/d) 2 
of the central star's intensity (g is the wire thickness). Note that (D/d) is the number of 
wire spacings into which the entrance pupil is sliced in each direction along the grid's axes. 

5.2. Observing strategy 

Our theoretical results suggests that fiducial spots should be placed in areas of the 
highest contrast of the image, where there are fewest remnant speckles, if the most accurate 
astrometry is desired. However, we recognize that they would then occupy a part of the 
optimal search region in the coronagraphic image. An appropriate observing strategy can 
remove this problem, or the observer can simply choose to ignore the region under the spots. 
If the spots are placed outside the AO control area and a significant amount of speckle 
noise is expected around the spots, thicker wires should be used in order to increase spot 
brightness. That way accurate astrometry could still be performed using the spots. The 
pupil mask can remain in place during observations, and yet not obscure the darker regions 
of the focal plane. Detectors that are destined for use in the next generation of ExAO 
systems typically utilize 2048-square arrays. If these detectors are used in purely imaging 
modes they can cover a wide enough field of view to enable the satellite spots to fall outside 
the AO control area. However, if an integral field spectrograph configuration is used, the 
number of pixels required to cover the dark region of the focal plane in both spectral as well 
as spatial dimensions becomes a concern with these detectors. In such cases the satellite 
spots can be placed within the AO control area. 

The astrometric device we describe could be left in the light path of the instrument all 
the time, or it could be removed from the beam, to be reinserted as needed. The presence 
of these fiducial spots may hamper the study of faint diffuse structure around the target 
star, and may also occupy too much of the best search areas in the coronagraphic image. 
The reticulate mask must be positioned at an angle that is not coincident with the telescope 
spider structure, or the calibrator spots may be lost in the glare of the diffraction spikes, 
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and the spots' morphology confounded by residual spider diffraction (Sivaramakrishnan & 
Lloyd 2005). 



The fiducials spots' locations and morphology will need calibration. Binary and triple 
star systems, or fortuitous groups of stars with appropriate brightnesses will need to be 
observed. Given that the light diffracted by the grid will diverge from the light of the parent 
star, starting at the location of the grid, a consideration of geometrical image distortion is 
merited on a case-by-case basis for each instrument that uses this device. There may be 
some utility to performing astrometric calibration of the spots without the coronagraphic 
focal plane mask in place, but with the undersized Lyot stop on classical coronagraphs. 



With four fiducial spots in the instrument's field of view, the spots provide most infor- 
mation about the star location if they are treated together. Due to the spectral smearing in 
the radial direction, each set of two spots lie on a line intersecting the central star. These 
two sets of two spots can be used independently to pinpoint the star in each of the two or- 
thogonal k x and k y directions, in the absence of atmospheric differential refraction. In actual 
data reduction, we suggest fitting a line to each set of two radially smeared spots. The fits 
of these spots should be achieved in a fashion similar to that used for optimal extraction of 
spectral data from digital images, (e.g., Oppenheimer et al. (1998), and references therein.) 

The intersection of the two lines from the two sets of fits pinpoints the star. This 
technique is independent of the spectrum of the star, since for each fit only the shape of the 
smear PSF in the direction perpendicular to the spectral smear is used. As a scientific bonus, 
a very low-resolution spectrum of the star can be recovered from each observation using the 
reticulate mask. The ultimate precision of this technique is limited by photon noise. 

The uncertainty in judging where the center of two radially-smeared spots lie (in the 
direction perpendicular to the spectral bandwidth-induced smearing) is of the order of the 
size of the pattern in this direction — viz., a spatial resolution element — divided by the 
square root of the number of photons in the two satellite images (assuming the spatial 
resolution element is at least sampled at the Nyquist frequency): 



Here the stellar spectrum is F£, and the detector has a quantum efficiency as a function of 
A of Q\. h is Planck's constant, and c is the speed of light in a vacuum. This assumes that 
the spots are significantly brighter than the surrounding speckles. 



5.3. Relative astrometry 




(12) 
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5.4. Relative photometry 

In conducting the optimal extraction fits of the smeared spots as described in the pre- 
vious section, a very coarse spectrum of the star is retrieved four times. The sum of the 
counts in these four spectra is 4e 2 multiplied by the stellar intensity (unocculted by the 
coronagraphic mask, but including the throughput effects of the coronagraph) . The derived 
intensity can then be used in all images with the reticulate mask in place to calibrate any 
other object in the field of view, either relative to the primary star, or in an absolute sense 
if the primary star is well calibrated and measured separately as a flux standard. 

6. Discussion 

We have shown that a robust, easily- made optic can be incorporated into existing coro- 
nagraphs, or designed into future ones, in order to create controlled fiducial spots in the 
coronagraphic image plane for quantitative calibration. This device is easy to use (in the- 
ory), and, furthermore, does not require changing the phase of a wavefront with a deformable 
mirror. The latter quality could simplify the way coronagraphic astrometry and photometry 
on a space-based optimized coronagraph is performed. At the limit, on the planned ExAO 
coronagraphs being constructed over the next five years or so, we expect about 1/100 of 
a pixel in astrometric accuracy for the primary, and 2% photometry. The accuracy of the 
relative astrometry between the apparent companion and the coronagraphically-suppressed 
central star depends entirely on the statistical significance of the companion detection. Fidu- 
cial spots placed inside the region of the image plane controlled by the deformable mirror 
are best for astrometry; the optic that generates the spots can be removed and replaced as 
needed since astrometric sensitivity to tilts and lateral shifts of the optic are low. Labo- 
ratory and on-sky use of such a device will develop our understanding of its behavior and 
utility. Marois et al. (2006) present early laboratory test results on the behavior of a similar 
type of fiducial spot. We plan to implement our reticulate pupil mask in the Lyot Project 
coronagraph in the near future, for on-sky tests of the utility of this technique. 

If the wires in the grid are of the order of a few wavelengths of light wide, these simula- 
tions may not apply, and one must abandon the Fourier optics treatment and use the Fresnel 
approximation or full electromagnetic propagation methods to calculate the expected loca- 
tion and characteristics of the fiducial spots. Curvature and distortion of optical planes 
- for instance, the image of the pupil where we place our grid — might require our grid 
geometry to depart from its simple periodicity in order to create compact, well- understood 
satellite spots. For extremely demanding applications in future instruments these issues will 
have to be addressed. 
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Fig. 1. — Monochromatic and broad-band direct and coronagraphic PSFs with a square- 
geometry reticulate pupil mask. All images are on a logarithmic greyscale stretching 10 
magnitudes fainter than their peaks. The pupil is 128 pixels across, and the grid has a wire 
spacing of 16 pixels, with 2-pixel wide wires. 1: direct PSF for the shortest wavelength of 
a 20% bandwidth filter with uniform transmission within the bandpass, in the absence in 
of phase errors. The satellite PSFs off the origin but along the horizontal and vertical axes 
are fainter than the central core of the PSF by a factor e 2 = (g/d) 2 , where g is the wire 
thickness and d the wire spacing. The satellite spots off the axes are e 4 fainter than the 
corresponding central peak. 2: coronagraphic PSF at the shortest wavelength of the filter. 
The off-axis sea of satellite spots are more visible in the coronagraphic image because the 
core has been suppressed. 3 and 4: direct and coronagraphic PSFs for the longest wavelength 
of the filter. 5 and 6: direct and coronagraphic PSF for the full bandpass. The length of 
any particular radial streak in this last pair of images, (in resolution elements at the central 
wavelength of the bandpass) is approximately the fractional filter bandwidth multiplied by 
the radial distance of the spot at band center. The streaks all point towards the origin, so 
the smearing has no effect on astrometric precision according to Fraunhofer regime image 
formation theory. We suggest using the four satellite peaks closest to the core as fiducials 
for the position of the central occulted star in coronagraphic images. 
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Fig. 2. — The essential planes and stops in a coronagraph. The entrance aperture is A, 
the direct image at B falls on a focal plane mask whose transmission function is m(k). 
The re-imaged pupil plane D, after being modified by passage through a Lyot stop with a 
transmission function N(x), is sent to the coronagraphic image at F. We place a square grid 
of opaque wires over the pupil plane A to create controlled fiducial spots in the coronagraphic 
image at F for astrometric and photometric purposes. 




Fig. 3. — Monochromatic satellite spots in the coronagraphic image of a band-limited Lyot 
coronagraph. The origin of an exactly periodic wire grid in the pupil was placed at 9 locations 
in a unit cell, chosen at random. Completely symmetric satellite spots result, regardless of 
the in-plane placement of the grid. The logarithmic scale is 5 astronomical magnitudes for 
the central panel, and 10 magnitudes for the surrounding panels. The grid was modelled 
with 16 wires, each two pixels wide, placed across a binary circular unobstructed 512-pixel 
diameter pupil. This wire is much thicker than that needed for a coronagraph producing a 
dynamic range of the order of 10 6 or more. 
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Fig. 4. — The effect of wavefront errors on astrometric fiducial spots. Monochromatic satel- 
lite spots in the coronagraphic image of a band-limited Lyot coronagraph are calculated 
for wavefront errors that would produce a Strehl ratio of 95%. In each of the three rows 
(labelled A, B and C) the origin of an exactly periodic wire grid in the pupil was placed at 3 
locations, chosen at random, in a unit cell. The "No WFE" column shows resulting corona- 
graphic PSFs on a logarithmic stretch between the peak satellite PSF intensity of 3 x 10~ 4 
(of the reference unocculted PSF using the the corongraphic Lyot stop), and a factor 251 
times fainter (6 magnitudes) below the satellite PSF intensity peak, and 14.8 magnitudes 
fainter than the central occulted star. The simulated AO system (with spatially-filtered 
wavefront sensing) used in the N act 24 column has 24 sensing and control channels across 
the entrance pupil diameter, and the N act 40 column images result from AO correction with 
40 channels across the pupil, with the phase rescaled to possess the same variance as the 
N act 24 case. The reticulate grid in the entrance pupil was modelled in the same manner 
as for Fig. 3. Residual speckles from uncorrected phase aberrations complicate the spots' 
morphology if the spots lie outside the dark square of the AO control region. This suggests 
using the right-hand arrangement, where fiducial spots are placed inside the AO control area 
for the most accurate astrometry. 



